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SUMMARY 


On the basis of a kinetic hierarchy system developed previously by 
the authors, fluid-dynamic moment equations are derived governing the inter- 
action between turbulent and thermal fluctuations. Truncation of the hier- 
archy under the hypothesis of "ternary" molecular chaos provides a unique 
closure condition on the set of equations. The (multipoint) kinetic theory 
is shown to reduce the inherent complexity of the conventional formalism 
of compressible turbulence theory and to minimize arbitrariness in formu- 
lating the closure condition. This assertion is based on two factors 
characteristic of kinetic theory: First, all the turbulent correlation 

terms appearing in ti. ' averaged equations of viscous compressible flow are 
expressed as expansion coefficients of a single function, viz., the two- 
particle correlation in the 2ji-spac3. Second, the closure conditxon is, 
then, equivalent to deriving and solving the equation governing the two- 
particle correlation. Actually, the first feature reflects favorably on the 
fact that turbulence corrections to fie gasdynamic equations representing 
the Reynolds stress and the turbulent heat flux, are each formed by single 
terms. Each of the expressions remain invariant in form whether the flow 
is compressible or not; such invariance is not to be expected in the conven- 
tional formalism of turbulent correlations. A superficial similarity of our 
formalism to a formalism based on the concept of mass-averaging breaks down 
at the second point at issue where our closure equations have a unique form 
without recourse to modelim. An S posteriori way of reproducing the whole 
set of equations (except for thermal-agitation terms) on a proper reinter- 
pretation of the conventional fluid-dynamic equations is indicated. In 
fact, a bilinear transformation of fluctuating quantities shows that the 
two formalisms are identical with regard to terms of double correlations. 

The essential difference between the two systems lies in the appearance of 
higher order correlations; turbulence corrections by the kinetic theory 
★ 
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include no such terms in principle, in contrast to the triple and the quad- 
ruple correlations intervening in conventional expressions of Reynolds 
stress and turbulent heat flux, respectively. 

1. INTRODUCTION 

1.1 Macroscale and Microscale Fluctuations 
in Nonequilibrium Gases 

Due to its enormous number of degrees of freedom in motion at the 
molecular level, a macroscopic body of a gas exhibits various kinds of 
fluctuations around each "mean" value, viz., each thermodynamic variable 
of the gas, such as the density, the pressure, the fluid-dynamic velocity, 
etc. Roughly speaking, we can subdivide the whole fluctuation into two 
classes; micorscopic and macroscopic fluctuations. The microscopic fluc- 
tuation, for example, the thermal agitation, is characterized by the fact 
that its relative intensity has the form of a(6N)”^/^ , where 6N is the 
number of molecules in the body of a gas under consideration, and a is 
a certain function of the thermodynamic quantities. This quantity is, then, 
vanishingly small in the continuum limit where thermodynamics and fluid 
mechanics are spoken of except under anomalous conditions (e.g., the entropy 
or the specific volume fluctuations at phase equilibrium; a —*■>») , On the 
other hand, a feature of the macroscopic fluctuation is that its intensity 
is independent of 6N. Therefore, in the continuum limit, this fluctuation 
alone is detectable by means of macroscopic measurements. Another feature 
of the macroscopic fluctuation is that its existence is limited to the non- 
equilibrium state, whereas thermal agitations prevail in the equilibrium 
state as well as the nonequilibrium state. Among wel^-known fluctuations 
of the macroscopic type are fluid-dynamic turbulence, unstable phenomena 
in plasmas, and the mild (low temperature) ignition of a hydrogen/oxygen 
system. 

In the present work we confine ourselves to a study of one of the 
interactions between microscopic and macroscopic fluctuations; namely, 
that between microscopic thermal fluctuations and the macroscopic fluctua- 
tions due to fluid-dynamic turbulence. Our objective is to derive equations 
governing the interaction of the basic turbulent field with the two fluc- 
tuations in the fluid-dynamic space for the case of a classical ideal gas 

* 

in first-order translational nonequilibrium. 
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1.2 Differences in Generation and Propagation 
of Fluctuation Correlations 

The microscale and macroscale fluctuations differ not only in their 
dependence on the scale of the system but also in the ways they generate 
and propagate correlations. The thermal fluctuation is generated by a 
spontaneous deviation from the mean value charaterizing the equilibrium 
state and is dispersed by relaxation. Microscopically, the generation and 
propagation of the fluctuation is via molecular co"? ision, and the corre- 
lation formed by this mechanism relaxes as^ exp(-c^t/e) , where € is a 
small parameter proportional to the mean free path, and c is the iso- 
thermal speed of sound. Then the characteristic time for the relaxation 
is tj ~ e/c^, a time proportional to the mean free path. 

Only recently , the existence of a second mode of propagating correla- 
tion has been revealed in addition to the one via direct molecular colli- 
sions: Correlation generated and propagated by collective motion of macro- 

2 3 

scopic bodies of the gas. Grad and Sastri state that in contrast to the 
direct correlation propagation, which persists only for a few mean free 
times, the second mode of correlation must propagate at least over a hydro- 
dynamic time, judging from the fact that the correlations are shown to 
obey certain "conservation" equations of two-point fluid mechanics. In 
fact, decay of the second mode of velocity-velocity correlation in a gas 
at rest obeys a diffusion equation and therefore behaves like exp(-x /et). 
This second mode, then, has the time constant t^^ ~ L^/e (where L is 
a hydrodynamic length) and therefore decays much more slowly than the first 
mode with t^ ~ e/c^. This implies a striking feature of the second mode: 

It must be detectable by macroscopic (fluid-dynamic) measurements. Further- 

4 

more, it has been shown that the correlation of the second mode is even 
amplified if the flows under consideration are unstable with respect to 
the cl?ssical criteria of hydrodynaraical stability. Fluctuation evolution 
in a blafius boundary layer, calculated along the line of two-point kinetic 
theory, has exhibited agreement with the experimental results of Schubauer- 
Klebanoff. Hence, we conclude that the second mode of correlation has a 
certain implication as a tuioulent correlation. 
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1.3 Description in Terms of 
Microscopic Density 

A crucial factor of a kinetic theory in dealing with fluctuating 
phenomena is the quality of statistical processing that has to be invoked 
to make the problem tractable. For example, it has been shown that instan- 

4 

taneous turbulent fluctuations are wiped out in defining the Boltzmann 
function as an expectation value for the number density. Also, the thermal 
fluctuations are missed^ even in the most general formulation of the BBGKY 
formalism. To avert this difficulty, it is advisable to start with an 
exact distribution function where no probability concepts have yet been 
invoked and where therefore every fluctuation is self-contained. To this 
end, we shall start with a distribution which has been proposed in plasma 

g 

kinetic theory, and has become known as the microscopic density. It is 
defined by 


N 

f (z,t) = m ^ 6 ^z - z (t)j 
s=l 


( 1 . 1 ) 


where m is the mass of a particle, and z = denotes a phase-space 

point, (t) gives the locus of the s^^ particle in the phase space, 

6 is Dirac's delta function, and the summation is over all particles N 
under consideration. The microscopic density has a favorable feature in 
comparison with the Liouville density function which has heretofore formed 
the basis for the BBGKY hierarchy: The microscopic density is defined in 

the (six-dimensional) phase space, in contrast to 6N space of the Liouville 
density, and has a definite physical meaning such that the expression 


be Av 




o 

f 


dx dv 


Ax Av 


gives the "exact" mass density at an instant t in the phase space z. 
This is an immediate consequence of the fact that each integrated delta 
function (1.1) is unity or zero de .iendir on whether the given particle is 
located inside or outside, respectively, the volume Z\x Av. 
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In terms of this distribution, macroscopic variables are defined such 


that 



(J-.2) 



I 


O —¥ 

V. f dv 

X 


m 


L 





(1.3) 


give, respectively, the density and the mass flux density which are subject 

to thermal as well as turbulent fluctuations. However, the presence of the 

s ) 

delta function and the unknowable nature of initial values of x and 

s ) o o 

V make it difficult to solve for p and m^ with instantaneous fluc- 
tuations restored (except by numerical experiments, e.g., the Monte Carlo 
method) . The simp] ast way to make the formulation tractable is to intro- 
duce a smoothing or av'.raging procedure whose operation on the microscopic 
density yields the .Joltzmann function, or more exactly, the one-particle 
distribution function of the hierarchy 



(1.4) 


In other words, f is the mean mass density in the phase space. Since, as 
noted above, neither turbulent ncr thermal fluctuations survive the average, 
all the instantaneous fluctuations must be included in 

Af = f - f (1.5) 


Accordingly, the density and the mass- flux density fluctuations are given, 
respectively, by 


Ap 



I 


Af d!v 


o 

Am. = m. - m, 
111 


^ v^ Af d"v 


( 1 . 6 ) 

(1.7) 
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where p and are the average density and the average rnass-flux density 

obeying the (Boltzmann) rormalism of the classical kinetic theory: 

( 1 . 8 ) 

( 1 . 9 ) 

Statistical behavior of the fluctuations is spoken of only through 
their correlations; for example, a tensor defined by 


p ■ j f dV^ 

. =• I V . f d”v 


Am. Am . 
1 3 




v. V , 

1 3 


Af Af d\T dv 


( 1 . 10 ) 


describes the average characteristics of the fluctuation in the mass-flux 
density at a seven-dimensional point (x,x,t) . In order to calculate the 
factor 


Af Af 


oo A 
ff - ff 


( 1 . 11 ) 


appearing in (1.10) , a simple mathematical identity 



I 

s=l 


needs to be noted. 
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Let the average, defined by fl.4), be taken over the identity; then along 
with definition (1.1) the following equation results; 


f(z)f(z) - fjj(z,z) +ra6(z - z)f(z) 


( 1 . 12 ) 


where f is defined by 

fjj(z,z) = 6[z - 2 ^®^ (t)] 6[z - (t)] (1.13) 

s^s 


is called the two- t>ar tide distribution function in the BBGKY hierarchy. 
Customarily, the ♦imction f^^ is decomposed further as 

fll(z,z) » f(z)f(z) + ^(z,z) (1.14) 


where f is called the tworparticle correlation. Then equation (1.11) 
reads, in view of (1.12) and (1.13), 


Af(z) Af(z) » + m6(z - z)f(z) 


(1.15) 


The physical meaning of each term on the right-hand side can be seen from 
(1.10) with (1.15). Assuming f to be a Maxwellian distribution, we have 




Af aJ d^ dv 


JviV^^(z,z) d'v dv 


+ 6 (x - jT) (pa) 


(1.16) 
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1 /2 

where n is the mean number density and a * (RT) is the isothermal 

speed of sound. Then the second term on the right-hand side of equation 

, 7 

(1.16) turns out to be the thermal agitation, whereas the first term 
represents the macroscopic fluctuation correlation in nonequilibrium situ- 
ations that is directly connected with the turbulent correlation.^ Thus, 
the macroscopic and microscopic fluctuations are uncoupled in (1.16), but 
they will be shown to interact with each other through the equations of 
fluid dynamics. 

In Section 2 we summarize our methodological basis for what follows; 
a hierarchy system which is based on the master Boltzmann equation instead 
of Liouville's master equation of the BBGKY hierarchy. Specifically, we 
discuss the one-point and two-point equations of the hierarchy which, to- 
gether with the assumption of "ternary" molecular chaos, form a closed set. 

In Section 3, factors causing a crucial difference from the conventional 
turbulence formalism are extracted and discussed. Fluid-dynamic moment 
equations out of the one-particle equation are described in Section 4 with 
special reference to the coexistence of compressibility and turbulence. 
Differences in several nontrivial aspects from the Chapman-Enskog ' s or Grad ' s 
classical schemes are pointed out. In Section 5, the connection of these 
equations with those derived on a macroscopic basis is discussed: The two 
formalisms are shown to yield the same equations in the lowest order 
correlations. In Section 6, the two-particle equation which plays a central 
role in the closure problem is detailed in its equivalent form of moment 
equations of two-point fluid dynamics. These equations are shown to be 
reconstructed except for thermal agitation terms on the basis of phenome- 
nological gasdynamic equations if dependent variables are identified with 
those defined in terms of microscopic density (Section 7). 

2. A HIERARCHY SYSTEM ON THE "MASTER" 

BOLTZMANN EQUATION 

2.1 Master Boltzmann Equation 

As is well-known, Liouville's equation, the master equation of the 
BBGKY hierarchy, has an implication as the equation of continuity in the 
(6N-dimensional) r space. A parallel procedure of generating a hierarchy 
system in the (6-dimensional) p. space turns out to be practicable, starting 
with the equation of continuity in this space. In reference 1 such an 
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equation for classical ideal gases has been shown to have the form of the 
Boltzman equation in which the Boltzmann function f is replaced by the 
microscopic density % defined in (1,1); 


o 

B(f) 



- 0 


( 2 . 1 ) 


where J is the conventional Boltzmann operator for the short-range xnter- 
molecular force 

a > 3 (2.2) 

and where z stands fcr a field particle on which the integrcition is 
effected. Equation (2.1) may be called the "master" Boltzmann equation as 
distinguished from its averaged version (eq. 2.3 below) governing the 
Boltzmann function. 

A salient feature of the master Boltzmann equation (2.1) is that no 
statistical processing such as the molecular chaos assumption nor even the 
probability concept (e.g., the distribution function as an expectation 
value for number density) has had to be invoked in deriving the equation.^ 

It is simply an exact continuity equation, in the p, space at each instant, 
of fluctuating number 




_lO 

^f dz 


of molecules in a specified volume 6z. The temporal resolvability of this 
equation is limited by the nature of the operator J that is unable to 
resolve temporal evoxution within the collision duration time t = r/a, 
where r is the characteristic radius of curvature of a field particle at 
molecular encounter. Therefore, t is interpreted as the time-resolvability 
limit of the equation. It is vanishingly small for a rigid spherical mole- 
cule (a -► “) , c.nd increases as a is decreased (softer collision) . We can 
show that the thermal agitations survive this fluctuation-filtering action 
of the operator J on condition that t be smaller than the mean collision 
time. In fact, spontaneous stress fluctuations caused by thermal motion 
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of a gas are shown to relax as^ exp(- pt/n) , (p , the pressure; (x, the vis- 
cosity coefficient) , so that the characteristic time t* of thermal agitw^- 
tions, T* ~ m/p ~ where e denotes the mean free path, is seen to 

be of the order of the mean collision time. Then the condition t* » t, 
or equivalently i » r (the necessary and sufficient condition for binary 
molecular collision) , ensures that thermal agitation will be rt •• *^ed by 
equation (2.1) for a certain class of "hard" molecular encounters including 
Maxwellian molecules (cr » 5) . 


2.2 A Hierarchy System 

In this section we present one-particle and two-particle kinetic 
ec[uations of a new hierarchy system on the basis of the master Boltzmann 
equation (2.1), with emphasis on the contrast with their BBGKY counterparts, 
(We do not duplicate here a derivation of the equations, which is carried 
out in some detail in reference 1.) The equation of the lowest order, viz., 
the one-particle equation, is obtained by averaging equation (2.1); 



(2.3) 


where (1.14) and (1.12) are incorporated, and where, on the right-hand 
side of (2.3), we have util^jed the fact that j|^6(z - z) fj =0, implying 
that no effects due to thermal agitations are observed in the averaged 
equation. In fact, this equation is exactly the same as the one-particle 
equation of the BBGKY hierarchy, and has the form of the Boltzmann equation 
with no molecular chaos assumption invoked. Evidently this equation differs 
from the classical Boltzmann equation by the presence of the term ^ in 
the collision integral so that the equation is not closed by itself. We 
need, therefore, to address ourselves to the next equarion of the hierarchy 
that governs the function This inherent difficulty of indeterminacy 

is not peculiar to the one-point equation. The two-particle equation 
includes higher-order terms as well; the three-particle distribution 
function fju intervenes in the equation in an unavoidable fashion. Then 
a question arises: What is the simplest possible form of truncation rule 

which retains the essential features of turbulence (•<// 0) ? Two points need 

3 

to be noted here. First, Sastri has pointed out that a truncation rule 
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( 2 . 4 ) 


be chosen arbitrarily, but should obey 


Aft::*'- 


to k 


J^Lij(fii,G)dz ■ CLj(f,fjj) 


(2.5) 


£'i^isit€ C is a constant and 


+ 1 ^ 


(N: I,II,III,,..) 


( 2 . 6 ) 


,th 


is j^the kinetic equation of N “ BBGKY hierarchy that is exact in the sense 
that idle Liouville equation is so. Condition (2.5) limits the choice of 
ta an euftent that the arbitrariness introduced in assumption (2.4) , 
thereby intervening in the two-particle equation does not "contaminate" the 
exactness of the one-particle kinetic equation. The second point at issue 
is that every G so far obtained and eligible in the sense above leads to 
an identical form when truncated at the three-particle level, whereas they 
differ at the four- (or more) particle level in effecting truncation. In 


fact, among the truncations meeting condition (2.5) are the well-known 
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cumulant truncation rule and the hypothesis of N-particle molecular chaos 


N 


% = n Af{z^) = 0 


(2.7) 


n=l 


where Af has been defined by (1.5) and where terms of macroscopic corre- 
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lation alone are taken into account. Both rules as well as some others 
are shown to yield an identical truncation rule at the three-particle stage 


f^jj(z,z,2) - fff - f^fz,2) - fi//{z,z) - fipiZyZ) = 0 (2.8) 


Thus, we have seen that the truncation at this stage and of this form is 
fairly general, although not perfectly so. 


Once the principle of the closure has been established, constructing 
an equation governing i}/ is rather straightforward. The two-particle 
equation is provided by 


Af(z)B[f(z)] + Af(z)b[f(z)] - 0 


(2.9) 


which, in view of (2.1), (2.3), and (2.8), is written in an alternative 
form^ 

^ V'(z,z) - p(z)^(z'*',£)j 


- K 12*^151 |f (2) V'(ztz)j = xj^z'^zj 1^6 


(z - z)fjj(z 




k[j+|s][, 


5(z - z)fj.j(z 


,z'*‘)J - 6(z - z)Jp'^zj jfji(z»2^)j (2.10) 


with the operator K defined by 


K(z*^jz)g(z,z'^,. . .) = J(z^jz) |^g(z,z'*’,. ..) + g (z'*’,z , . . . )j (2.11) 


where z^ indicates a field particle on which int ^gration operates. If 
the terms on the right-hand side of equation (2,10) are set to zero, the 
equation reduces to the two-particle equation of the BBGKY hierarchy that 
is homogeneous in In the case of a non-zero right-hand side, note that 

the inhomogeneous terms have a common factor 6 (x - x) , which implies that 
they stand for effects due to thermal agitations. It should be stressed 
that equation (2.10) is exact witli regard to the terms arising from thermal 
agitations; The closure approximation has been incorporated only in the 
macroscopic part of the three-particle distribution function (see eq. 2.8). 

Equations (2.3) and (2.10) constitute a closed set of governing equa- 
tions for f and ^ where the effects due to thermal agitations, missing 
in the BBGKY formalism, are taken into consideration. The remaining part 
of the present work is devoted to a derivation of the fluid-dynamic moment 
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equations which are correct to the first-order in (double correlation) , 

the mean free path e, viz., the order of magnitude of transport coeffi- 
cients, and n ^ {n: the mean number density) representing the first-order 

effect of thermal agitations. 

3. EFFECTS OF TURBUI£NCE IN THE 
KINETIC THEORY FORMALISM 


3.1 Expansion in Hermite 
PolyncHnials 

In order to accomplish transition from. the kinetic to the fluid-dynamic 
equations we expand f and in a single and a double series of Hermite 
polynomial respectively, and discuss equations governing their expan- 
sion coefficients. A systematic method of moment expansion has its origin 

9 

in the so-called 13-moment method of Grad. Since the expansion is around 
the (local) Maxwellian distribution (eq. 3.3 below) expansion in three- 

dimensional Hermite polynomials is considered to be most relevant. The 13- 

9 4 

moment expansxon for f and its two-particle counterpart for f are the 

following: 

f(z) = (3.1) 


^(z,z) 


= (pp)"'f 


^(z)fo(S) 


(3,3) 


(J,K) 


(J,K) = (0,0) 


J*K 

c c j;r,: 


j . . k. . 


(3.2) 


with 


f^iz) 


p (2vc^) 




(3.3) 


^ = c~^ (v - v>) 


(3.4) 


where ^ is the Hermite polynomial defined in terras of the non-dimensional 
peculiar velocity also p, p, c and "w are the stress deviator 

tensor, the heat flux density, the density, the pressure, and the two ref- 
erence velocities to be determined later. Each of these quantities except 
p is to be generalized properly to include effects due to turbulence. 


13 


Two remarks in connection with expression (3.1) are in order at this point 
First, we postulate the turbulent mean pressure p to ha^/e the same mathe- 
matical expression as it does in the laminar case regardless of its physical 
relevancy , 

P “ j + 3^ f d7 (3,5) 

/ 2 \ ( 2 \ 

where * - 3 is the contracted Hermite polynomial of the 

second order. Second, we postulate also that 

p^j = 0 (3.6) 

be retained together with the condition of vanishinc first-order terms . 

(See 3.L) Postulate (3.5) in which (3.1) is substituted for f deter- 
mines the parameter c as 


c 


2 


P 

P 


(3.7) 


Also, expansion (3,1) substituted in (1.9) determines the parameter w of 
(3.4) as 


w. 

1 



(3.8) 


With regard to expansion (3.2), there are no such degrees of freedom re- 
maining, so that the summation in (3.2) should span over all possible terms 
including nonvanishing first-order as well as contracted second-order terms. 

All of the expansion coefficients in (3.1) and (3.2) are solved in view 
of the orthonormality property of the Hermite polynomials as 


D. . 


‘ID 


= f dv 

J 


(3.9) 
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(3.10) 


!i 


3 ) 


f dv 


where 


^ijj 


5) 


and 


(X,?) = dv d? 

j • • )K •• I D** K»» 


(3.11) 


3.2 Instantaneous Fluctuations 

Expressions (1.8), (1.9), (3.5), (3.9) and (3.8) for p, , p, P^j> 
exhaust the set of 13 mean variables necessary for describing a fluid 
state. To each of these variables corresponds an instantaneous fluctuation 
some of which have been already given by (1.6) and (1.7). The whole set of 
the instantaneous fluctuations is given by 


% 


r- 

1 

Am. 

1 


w . + c 
1 1 

Ap 

r 

-J 

c- 

APij 


IT 

AQ. 


2 X 

L 




Af dv 


(3.12) 


whose form should be self-explanatory in view of the corresponding expres- 
sions for the mean values. 

It should be noted that all the variables appearing in equation (3.12) 
have in common the characteristic of being proportional to number density. 

In fluid mechanics, however, there is another class of quantities whose 
m.embers are independent of number density, for example, the velocity or the 
temperature. We have learned that quantities proportional to the density 
are more fundamental than those which are density-independent in the sense 
that fluctuations of the former quantities are directly defined in terms of Af. 
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It is sometimes required, however, to express the fluctuation formulas 
with respect to the density-independent quantities. For example, laser- 
doppler velocimetry picks up the velocity fluctuation rather tlian that of 
the mass-flux density. Similarly, some measuring devices perceive the tem- 
perature fluctuation instead of the pressure fluctuation. 

The density-independent quantities can be defined only indirectly. 

We shall show this for the velocity and the temperature. The instantaneous 
o 

velocity u. is defined as 


o 
u , 



(3.13) 


Separating this expression into average and fluctuating parts, we have 


m. = pu. + Ap Au. (3.14) 

Am^ = p Au^ + u^ Ap (3.15) 

As for the temperature, it would be most natural to define this quantity as 
equi-parti tioned kinetic energy per molecule due to purely thermal motion, 
namely, that excluding the turbulent energy. Then we have 



o — 

u . ) - f dv 
1 


(3.16) 


where k is Boltzmann's constant and n is the fluctuating number density. 
Splitting this expression into average and fluctuating parts leads; to 


pRT + R Ap AT = p 



(3.17) 


RT Ap + Rp \T = Ap 


(3.1B> 
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where 



(R 



and where formulas (3.5) and (3.12) have been used. 

3.3 Physical Meanings of 

J « * X, • • 

With these preliminaries on the instantaneous fluctuations, we can 

(J k) 

show a more direct physical meaning for R' ’ ' thafi is implied in (3.11). 
As is seen from (1.15), the R's are linearly related to 


AA = 


f 


a(z)b(z)Af Af dv dv 


(3.19) 


namely, the fluctuation correlation between two thermodynamic quantities 
A and B at physical- space points x and x. Formula (3.19) is an 
immediate consequence of averaged prodvict of <;-"pression (3.12) at the dif- 

A , A 

ferent points z and z. Quantities lV\ and AB stand for fluctuations 
of macroscopic quantities standing on the left-hand side of (3.12) , whereas 
a and b stand for those counterparts appearing within parentheses on the 
right-hand side. A straight-forward calculation of (3.19) in v/hich Af Af 
is replaced with the right-hand side of (1.15) yields 


Ap Ap = R^°=°^ + rnp 6 ("x - St) 


■A 


Ap 


+ p + mm^5(}T - 5) 


- 4 - m 4 . 


Am. Am, = -tm.p ‘R'A'* 

1 £ 1 , i 1^1 


+ p R . 

i 1 


( 1 , 0 ) 


(3.20) 


+ (pp) ^m. fn ,R + m(p ^m.- + p6 . ) i') (^ - x) 

1 lx 


Ap Ap = 1 + mpiMjT - ^) 


3 rr 




(Continued on next page ) 
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and 


1 -U.2) +i„-'„,R(0.2) 


iim.Z'p = -T R.' ' ' + y 
1 ^ J i,rr 3 


rr 


+ c^'r 


( 1 , 0 ) 


, 2 

+ c p m^R 


( 0 , 0 ) 


+ (x - x) 


The above formulas show that fluctuation correlations, in general, are com- 
posed of two factors: thermodynamic fluctuations, which are characterized 

by the presence of the delta function, and nonthermal, nonequilibrium fluc- 
tuations R, that are connected with a nonchaotic situation / 0) 

(c.f. 1.16). On the other hand, terms on the left-hand side of (3.20) are 
rewritten in terras of Au^ and AT by the use of (3.15) and (3.18). Term- 
by-term comparison of these expressions with those on the right-hand side of 
(3.20), where the thermal agitation terms are neglected as of higher order, 
leads to the following relationships; 


(0,0) 

<CL 

<1 

Q. 

<1 

il 

(1,0) 

i 

- p Au^Ap 

(1,1) 
i , £ 

= pp Au^Au^ 

(2,0) 

= 3pR AT Ap 

(1,2) 

i 

= 3ppR Au^.\T 

(2,2) 

= 9ppR"AT AT 




) 


(3.21) 




3.4 Notes on the Pressure Concept in 
the Presence of Turbulence 

We have introduced the average pressure p in the form of (3.5) for 
the sake of conformity with the expansion (3.1) which refers to the pecu- 
liar velocity f defined by (3.4). In fact, from the formal point of view, 
this seems to be the most natural definition of the pressure extended so as 
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to allow for turbulence. According to th la definition » however » p includea 
- effects due to ttttbulenc.e_a_s_well^^^^a^^^^^ motion of gas molecules, 

p f p 

o 

provided that p is to represent the instantaneous pressure iR:luding 
thermal effects alone, in order to obtain the relationship between p and 

S 

ir - p (3.22) 

we proceed as follows. A key to the correct form of p is provided by 
(3.18) which is exactly the fluctuating part of the ideal gas relation 

o oo oo 

p - RpT = knT (3.23) 


Since we have the definition (3.16) of the fluctuating temperature at hand, 
equation (3.23) leads to the definition of p as 



(3.24) 


Averaging this expression and utilizing equations (3.22), (3.21), (3.4), 
and (3.5), we have 


p = 7T + 


ifRd,!) 

3p i,i 

i- 


X 


(3.25) 


This relationship tells us that the quantity p defined by (3.5) includes 
effects due to turbulence in addition to the "genuine" pressure ir which 
is defined by (3,22) and includes thermal effects alone. With this rela- 
tionship, equation (3,17) is interpreted favorably: In fact, the equation 

reads 
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or 


TT = pRT + R Ap AT 



(3.26) 


which is nothing but equation (3.23) averaged. The difference between 
TT and p, which is ascribed to random motion due to turbulence, is of 
0(M®A®) in comparison with the pressure (where M is the characteristic 
Mach number). Therefore, the difference is negligibly small for incompres- 
sible flows, but can reach an appreciable magnitude in certain hypersonic 
situations . 

Correction formulas (3.25) and (3.14) tell us that the local Maxwellian, 
when written in terms of the average velocity and temperature, takes the 
form 


f 

o 


exp - 


V. - u - 
3 -J ^ 3 


RT + 

3p2 3p- 


All the correction factors appearing in the expression disappear for incom- 
pressible flows because of non- fluctuating density = o| and of 

vanishing Mach number R^^’^Vp^RT O(M^) « 1. 

3 j 3 

4. TURBULENT GASDYNAMICS ON THE BASIS OF 
THE CORRECTED BOLTZMAIIN EQUATION 

4.1 Moment Equations 

The only difference between the modified Boltzmann equation (2.3) and 
the conventional one is the presence of the term f in the colli -n 
integral. Thus, no difference is observed in the three moment equations 
(of conservation) corresponding to three summational invariants as moment 
functions. This situation is readily seen from the following expression 
of tha general moment equation for equation (2.3) with the moment function 
a(z) : 
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^j"af dv+j^Jav^f dv d7 . ij" [[a]] [f( 2 )f(z) 


+ y(z,z)j dK d"v dT (4.1) 


where the notation 


[[a]j - a(z') + a(z') - a(z) - a(z) (4.2) 

has been introduced. In deriving expression (“^.l) we have assumed that 
quantities appearing in the collision integral cf equation (2.3) are sym- 
metrical with respect to z and &. This condition is met with ip since 
we have put jT = x to secure collision encounter. The three conservation 
equations corresponding to the following choices 


a - 1, 


V 


2 


for which we have 


[[“E 

are shown to have the form: 

D 


^ (p'^m.mj + p6. j + p. .) = 0 


(4.4) 


3t 


f p + T p ‘"j 


dx . 
3 


ip‘Vm? + |pp’'mj + p"'ni.p. . + Q. 


= 0 (4.5) 


These equations are formally identical to their classical counterparts so 
far as the actual forms of and have been left unspecified. 

As has been seen above, the way in which nonvanishing ^ affects the 
fluid dynamic equations (4.3) through (4.5) is via p^^ and , whose 
actual expressions are to be determined from moment equations with moment 
functions ^nd , respectively. Since these moment functions 

are not sununational invariants of collision, the moment equations neces- 
sarily include effects due to V'* After some manipulation, we have moment 
equations, namely, equation (4.1) for a = and a = , 

respectively: 


^Pii , a 


m 


ao. ao. 

0 OQr-' 

■N 


f 



2 A ^ 

■7‘ij ^ 

Pir ^ 

P 


a 

mi 

2 A ^ 



m. 

1 



Pjr ^ 

P 

s 

p 


P 

dx. 

1 

P 


1 ^ 


^(p. . + i 

III riD p 


j^(2,0) _ j^(l,l) 

ij i I j 


4 6..r( 2,0) ^ 1 ^ j,(l,l) 

3 1 ] 3 1 ] r,r 


(4.6) 


and 


. a 
"sr 35 T 

r 


m 

.7 a 

m . 
1 

.2 a 

"*r 

.2 a 

m 

r 

— Q- 
P 1 

r 

P 

+ 5 — 

1 

P 

5 'll -s^ 

P 


+ RT 




ir 


. 7 ^ a 
357 * 7 Pir^ 


£i 

P 


^^(P6 +P)+4P-^ 

dx sr sr 2 ‘ jx. 


. i£B 



,(3,0) _ 2Rf2,l) + Rf2,l) 


ir ,r 


(4.7) 
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where B is related to the viscosity coefficient 
tivity coefficient > as 


4 


and the heat conduc- 


mRT 
^ “ 6B 


and 


(4.8) 


. 5mR^T 
^ “ 8B 


All K's in equations (4.6) and (4.7) represent the correlations at the 
same point (x = x) ; therefore, the relation holds. 

4.2 Generalized Navier-Stokes 
and Fourier Laws 

Equations (4.6) and (4.7) have elucidated, at the level of the 13- 
rnoment approximation, the manner in which the stress and heat fliix are 
affected by the presence of i}/, namely, the two-particle correlation. 
Neglecting the quantities of the order of transport coefficients in the 
left-hand sides of equations (4.6) and (4.7), that is, reducing the level 
of description to that of fluid dynamics, we obtain the following expres- 
sions for the stress and heat flux: 



^^ij^NS 



r( 2,0) +i6..R(2,0) 1 , ^(1,1) 

13 1,3 3 1 ] 3 1 ] r,r 


(4.9) 


and 


<°i>F 



-r!^ + 2r!^ 


p. 


( 2 , 1 ) 


ir,r 


(4.10) 


where ^^i^F given by the classical Navier-Stokes 

Fourier Jaws, respectively, as 


^^ij^NS 



and 


(4.11) 




Xi [pR 


Further simplification of expressions (4,10) and (4.11) is feasible if 

4 

we make use of the fact that the correlations between quantities containing 
a nonsummational invariant are smaller than those between summational invar- 
iants by a factor proportional to transport coefficients €. We have, for 
example , 


=16..r(2.0) 

1 ] 3 13 13 


r(2,0) 

_Li 

o(2,0) 


~ 0 ( 6 ) 


With this decomposition rule and the simplifying approximation, equations 
(4.10) and (4.11) read 


,p. ,) +1 rU.d .1 

p 1,] 3 13 r,r 


(4.12) 


- <° 1 >F 


(4.13) 


Substituting equation (4.12) into equation (4.4) with the relationship 
(3.25) between •rr and p in mind, we readily see that the term 
behaves just like the Reynolds stress appearing in the equation of incom- 
pressible turbulent flows. The "turbulent" Navier-Stokes relation (4.12) 
might look deceptively simple, however, when we are reminded that the rela- 
tion holds as well for compressible flows. In fact, we have imposed no 
restriction on the compressibility in the course of the derivation. Also, 
the fact that the turbulent heat transfer (the second term of 4.13) is 
expressed by a single term is not to bo expected in the conventional 
formalism. Actually, in the conventional approach, the averaged Wavier- 




m 


Stokes and the energy equations, in which mean and fluctuating parts are 
separated, are written in the forms: 


{pu.u. +TT6.. + p..^) =0 
dt 1 ox. 1 ] 1 ] ^ 1 ] 


^ [2 P“i^ + ph - t] + ^ |pu^ .“5“ 


with Pj^j^ 2^^ given by 


p. . = (p .)„c P Au.Au. + d.Ap Au . + u.Ap Au. + Ap Au . Au . 

ij NS ^ 1 j 1 ^ 1 ^ 1 j 


:4.14) 


(Q^)p + J pu^(Au^)^ - Y Ap Au^ f u^Ah Ap + p Au^Ah 


1 1 

+ h Ap Au^ + Ap Ah Au^ ~2 ^ “Au^ + -j pu^Ap (..'A’j) 


+ j Ap Au^ (AUj) 


(4.15) 


where h xs the specific enthalpy. These expressions for p. . and 

10 11 

are the axact versions of those appearing in tne textbooks on turbulence. ’ 


Note that the Reynolds stress includes triple correlation terms and that 
the turbulent heat transfer includes terms up to quadruple correlation. 

It is natural, then, to a.sk whether the two pairs of transport relations 
(4.12), (4.13), (4.14), and (4.15) are essentially or only seemingly 
different. This question is investigated in some detail in the following 


sections . 
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5. RELATION BETWEEN KINETIC AND MACROSCOPIC 
FORMALISMS OF TURBULENCE EQUATIONS 




5.1 Comparison of the Two Formalisms 

The system of ecjuations derived above, viz., equations (4.3) through 

(4.5) supplemented by (4.12) and (4.13), is. shown to have an intuitive 

interpretation in terms of S priori gasdvnamic equations whose dependent 

■ o o o o o 

variables are the Klimontovich variables, p, ni^/p) , p, defined in 

previous sections by 


o 
dm • 

L E + — 1 = 0 
O ot X . 

3 


(5.1) 




"j ~ "St dx^ 


o o 

o o 

-fei + g 6 + (P 

L p , 


= 0 


(5.2) 


L. = 


dt 


r o 21 


~ 002 OO « 

1 . V 

y - 1 P ^ 2 O 

L p J 

+ 

^2 02 -y - 1 O ^ 


o ,o 


+ + (8 




= 0 


(5.3) 


In the above, v is the adiabatic index and takes the value 5/3 for mona- 

o o 

tomic gases now under consideration, and F defined, 

respectively, by 


^Pjk^NS 


-M- 


^dx. - dx, - 3 


(5.4) 


(Q^)f = 


(5.5) 


We shall show that our kinetic turbulence equations are equivalent to the 
averaged version of equal’ ans (5.1' through (5.3) if, in the course of 
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( 


averaging, triple or higher-order correlations are disregarded and the 
transport coefficients (l. and >,) which depend on the temperature are inter- 
preted as material constants , and are therefore not subject to fluctuations, 
No explanation is needed for the continuity equation: Simple aver- 

aging of (5.1) gives (4.3). With regard to the Navier-Stokes equation, 
we proceed as follows. By the turbulent Navier-Stokes (TNS) equation wc 
mean hereafter momentum-conservation equation (4.4) in which p^^ is 
replaced with a.e generalized Navier-Stokes relation (4.12), namely, 


dm. 

— ^ -f- 
6t dx . 

3 


Tti.m • 
J- 3 


+ 'Pij’NS + p 




1 ,3 


=•■ 0 


(5.6) 


First, we note chat the limitation as to the equivalence of (5.6) and (5.2) 

o 

is not existent when the flow is incompressible (p = p: a material 

constant), in which case averaged equation (5.2) leads exactly to equation 

(5.6) in view of relation (3.14) with Ap = 0. For compressible flows, 

however, the identity is reached only conditionally because of the density 

variation. Then a straightforward averaging of (5.2), written in terms of 

m. instead of u. yields 
3 D 


i.h.s. of (5.6) + 0(e = 0 


(5.7) 


where A stands for fluctuating quantities. It is easily checked that 
the terms A^ and include a common factor Ap , so that equation (5.6) 

is exactly retrieved for incompressible flow. For compressible flows, 
identity of the two formalisms is fulfilled only within the approximation 


€ A^ = 0 
A^ = 0 


] 


(5.8) 


If, in the course of deriving (5.7) , the variable u^ is used instead of 
m^/p, we obtain equation (4.14). Apparently the form of (5.6) is more 
convenieiit than that of (4.14) in its invariance with respect to 
compressibility . 
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.Ji . ' 




- 4 

3 13 r,r 




(3.9) 


The former equation, with average taken and relationships (3.20) incorpo- 
rated, reads under approximation (5.8), 


— r p + m? + 
y - I 2p 3 dx. 

w J ^ 


^ m? + ^ m .mf + — ^ Y "c vP- • ^ mc 

2p 3 J 2 d“ ^ ^ Y - 1 P 3 p ] ^3 NS 


'! 


The same method is shown to be applicable to correlate the £ p riori 
energy equation with the energy equation of the present formalism. The 
latter equation, which is equation (4.5) with turbulent Navier-Stokes tind 
..•ourier laws (4.12) and (4.13) incorporated, is written as 

^ [3 „ . i ^3 ' 1 . _d_ V .iLd,!) 

^[2 P + 2 p J ^ax~[2p^ 2 P p + p y^j^NS +p [^,3 


p 1,3 3 13 r,r 


+ (Q.)„ T ^ — ~Y ~ 

'^^ 3 ' F 3 (7 - 1) p 3 


i r(I-.I) 

P 1,1' 


+ ^ 


4 

)- r .r 


1 r(1 ,2) 

P 3 


= 0 (5.10) 


where, in deriving the equation, t = p has been eliminated through the use 
of (3.25). This equation tells us that for a monatomic gas (y = 5/3) the 
equation reduces to equation (5.9) and further, that the factor 5/6 
appearing in the generalized Fourier law (4.13) is to be interpreted as 
"v/ 3 (y - 1 ) rn the case where the gas is not monatomic. 

It should be remarked that equation (5.9) (and eq. 5.10 as wall) is 
written in terms of the "pseudo" pressure p, so that it is net in perfect 
conformity with the turbulent Navier-Stokes equation (5.6) that is described 
in terms of the "genuine" pressure " and has a form easier to understand 
in connection with laminar cases. The parallel expression for the energy 
is feasible by introducing the internal energy per unit of volume: 
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E 





R.' 

1 


( 1 , 1 ) 


(5.11) 


Then equation 


(5.9) takes the form 


3 1 
2 



+ E 


N (m?m. 

\^2p" 


m • Vn. 

+ — *- (E + tt) + 

P F 


'Pij>NS 


+ i Rfiri'l 

P J 


+ 


(Qj)r 




0 (5.12) 


This expression is apparently consistent with that of (5.6) and also with 
the laminar case: It consists of physically identifiable terms provided 

that we interpret E + t as the enthalpy per unit of volume. This obser- 
vation allows the following interpretation for the turbulent gas. Turbu- 
lent flow of an ideal gas, in its average description, can be considered 
as an equivalent imperfect gas with thermal and caloric equations of state 
given by (3.26) and (5.11), respectively. 


5.2 Essen cial Difference Between 
the Two Formalisms 


As we have seen above, the kinetic and macroscopic methods have led 
to the same ecjuations of turbulence so far as the lowest order in corre- 
lation terms is concerned. However, it seems to be equally important to 
stress the difference between the two formalisms with regard to the appear- 
ance of higher-order correlations, because this point proves to be the 
essential feature of the kinetic- theory formalism that is not to bo attained 
by the conventional formalism. 


First of all, let us recall that the classical turbulent corrections 

to transport relations (4.14) and (4.15), which are exact within their 

regimes, include triple and quadruple correlations, respectively. These 

oo o oo O.- 

have their origin in the terms and in tl-.c Navier-Stokes 

and the energy equations, respectively. On the other hand, TNS equation 
(5.6) and turbulent energy equation (5.9), both of which are consequences 
of the kinetic theory, include only double correlations and 

R.' ’ ' through generalized Navier-Stokes and Fourier laws (4.12) and (4. 


through genei 

X 

respectively. It can be shown easily that the one-particle equations, 


13 ) , 
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namely, the Navier-Stokes and the energy equations of fluid mechanics, 
should not include correlations higher than double correlations, in 
principle. In other words, this is true independently of the multi-par- 
ticle stage at which truncation is effected. To prove this, only two facts 
need be remarked. First, a general property of the BBGKY hierarchy is 
that each hierarchy equation concerns only two distribution functions of 
consecutive hierarchies. In fact, the equation of each hierarchy has the 
form 


Df 

n 


(5.13) 


where D and J are certain streaming and collision operators, 
respectively. For n = 1, the equation reduces to equation (2.3). The 
other fact which is necessary is that triple correlation is directly 

connected with a nonreducible part of the three-particle distribution 
function 





in such a manner that 

.\A AB AC = J^a (z)h{z) c (z) dO dv 

therefore, nonvanishing terms of A^ are directly connected with nonvan- 
ishing appear in equation (5.13) for n = 1. Thus, we 

conclude, as a general rule, that terms of order A~^ do not appear in any 
of the moment equations stemming from equation (2.3). Actually, formulas 

(4.12) and (4.13) do not include such terms although no appro.ximation to 
rule out higher-order correlations has been effected. These retain a level 
of accuracy comparable with formulas (4.14) and (4.15) of the conventional 
expressions where triple and quadruple correlations intervene in unavoidable 
fashions . 

The above statement does not imply, however, that formulas (4.12) and 

(4.13) are not at all affected by the form of the function jjj (assumed 
or solved). It does affect them imp] ici tly through the double correlations 
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and whose actual forms are to be solved from equation (5.13) 

i 1 3 ^ 

with n = 2, where the term "'■''jjj appears on the right-hand side of the 
equation . 

6. CLOSURE EQUATIONS 

6.1 Comparison with Conventional Closure 
12 

M. Rubesin first pointed out the close similarity between the system 
of equations (5.6) and (5.12) and a system involving mass averaging, 

Mass averaging is an artifice which formally eliminates density fluctuation 
from the momentum and the energy equations, thereby enabling application 
of the modeling conventionally used in incompressible turbulent flows. In 
the present formalism an equivalent artifice is forbidden by the structure 
of the theory. Here, all the fluctuation terms appearing in the one-par- 
ticle system emerge from a single function as coefficients of its 

expansion, and the governing equation for ^ is already at hand (eq. 2.10)’ 
Therefore, invoking a special kind of averaging instead of the straight- 
forward one we have employed would make equations (3.25) and (3.14) over- 
determined. For the same reason, it would be redundant to introduce 
special heuristic arguments concerning the Reynolds stress of 

(4.12) and the turbulent heat flux of (4.13), such as, for example, 

the mixing length concept and a turbulent Prandtl number of unity. The 
evolution of these quantities hinges totally upon solution of the partial 
differential equations in (x,x,t) space derived from kinetic equation 
(2.10) by the moment expansion. 

In the case of incompressible flow, these equations have been worked 
out in reference 4. There, the equations are shown to be separable into 
tv;o groups of variables (>T,t) and (S,t), and thereafter reducible to two 
Orr-Sommer feld equations with different physical implications. The gene- 
ralization now at issue is two-fold; inclusion of compressibility winch 
allows for a nonvanishing and inclusion of thermal- fluctuation 

effects which give rise to inhomogeneous terms in the correlation equations 
The two-particle roment equations including these generalizations are 
developed in what follows. 


6.2 Two-Part.f cle Moment Equations 


A general expression for the moment equations of equation (2.10) has 
been derived by the authors (eq. (6,^) of ref. 1). The kinetic equation, 
after being multiplied by a moment function a(z){3(z) and integrated with 
respect to "v and v, is written in the form 









tii, 





vv 




with 







- j^|a(z) 3 (z)jj l^aCz) + a(z)j 

[p(z) + 3(Z)] [[a(z] 


fll(z,z) 



( 6 . 2 ) 


vv 


where the following abbreviations are employed: 


z d7, Z dK(z,z)dv* dv, 




\ (o.3) 


) 


£2 . ^ ^(Vp)Z 

£>t dt ax^ ’ Dt Ot 'x^ J 


It may be relevant to point out here that term G on the right-hand side 
of equation (6.1) gives the effects due to thermal agitations and that the 
term is nonvanishing unless the product a:- is a summational invariant as 
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well as <1 and P themselves. For example, the two-point fluid-dynamic 
equation describing momentum- momentum correlation is subject to an inhomo- 
geneous term ascribed to thermal agitations. These situations are actually 
the case in the following set of equations which exhausts the whole range 
of possib. choices of a and p when both are summational invariants; 


^ ( 0 , 0 ) 

H — = 0 
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P.Ei’RODUCiBlLITY UP' THE 
■miUINAL PAGE IS F(X)R 
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( 6 . 8 ) 
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with 


,(J,K) _ p(J,K) 


= R ' (x,x) 


(6.10) 
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In the above equations we have also employed the following expressions: 


(2,K) , r( 2,K) _ 1 g j^( 2,K) 
ij / ij / 3 ij 


^(3,K) _ p(3,K) 
r i _ R i 


( 6 . 11 ) 


( 6 . 12 ) 


where the moment ' corresponding to the second variable is a summa- 

tional invariant. Equations (6.4) through (6.9) are regarded as the evo- 
lution equations for correlations r' ’ ' between summational invariants. 
As the actual expressions show, these are not closed by themselves, but 
include correlation terms involving the r's associated with nonsummational 
invariants. In deriving terms on the right-hand side, viz., terms of ther- 
mal agitations, only linear terms in R or e have been retained. 

In order to have a closed set of equations we need to construct equa- 
tions for the r's; that is, the fluctuation correlations between the stress 
tensor/heat flux vector and summational invariants. To extract the essen- 
tial feature of equations governing the r's, an inductive approach seems 
to be more palatable than a deductive one. For this purpose we derive 
equations governing two correlations of (4.11) and (4.12) for K = 0: 
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These equations are simplified by imposing an order-of-magnitude estimate 


r ' 0(R‘,Rc) 


(6.15) 


which will be proved 5 posteriori. In fact, retaining only terms of 
0(1, R/e) and leading terms of delta function in the equations, we have 
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These are the two-particle versions of the Navier-Stokes and the Fourier 

laws generalized to include effects due to thermal agitations. These 

results also assure that postulate (6.15) actually holds. It should be 

remarked, however, that not all terms in equations (6.16) and (6.17) are 

equally important. For example, when (6.16) is substituted for 

equation (6.5) , all the terms except the first on the right-hand side of 

(6.16) find respective terms that are larger by the factor 0(e ^) or 

0(R~^). On the other hand, there are no terms to compare with regard to 

the first term of (6.16), therefore, this term alone should be retained. 

This is also the case with equa<"ion (6.17) , as well as with the expressions 
(2K) (3K) 

for general ^ ^ i* ’ ^^ving carried out the simplification, 

we are led to the following equations: 
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n 

The connection of these expressions with the Navier-Stokes and the Fourier 

laws of transport process should be self-explanatory. '■ 

6.3 Interaction Equations 

Equations (6.18) and (6.19) express fluid-dynamic correlations between 

noiisuminational quantities in terms of those between summational invariants; 

therefore, together with equations (6.4) through (6.9), they form a closed 

set. For incompressible flows (p, a materia] constant), terms including 

the density fluctuation disappear = oj , and the system reduces to 

a closed set of equations (6.5) , (6.5) , and (6.8) for K = 1. Eliminating , 

(12) 

r ’ from the latter equations, we have equations governing the inter- 
action of turbulence with thermal agitations in incompressible flows; 

= 0 (6.20) 
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where u. is the mean velocity and is identified with , *in. in this 
1 1 

situation. Equations (6.20) and (6.21) differ from the correlation equa- 

4 

tions of. the previous theory by the presence of the term with the delta 
function on the right-hand side of (6.21) This term, lacking in the BBGKY 
formalism, provides the only mechanism for generating correlation in a 
perfectly quiescent shear flow. Also, this term, which is ascribed to the 
effect of thermal agitations, is seen to be very small because of the 
presence of the factor n , which is of 0(10~ ) under normal temperaturt 

and pressure condition. Under realistic conditions, therefore, this effect 
in initiating turbulence is masked by other macroscopic mechanisms, such 
as noise or free-streami turbulence. Taking the thermal motion into account 
makes sense only in evaluating the maximum transition i.^ynolds number ;.n a 
shear flow. This evaluation has been attempted for the Brasius flow only 
qualitatively in reference 4, because the form of thermal agitation effects 
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Tr;P 


4 . / 


s not known at that stage. In this regard, eouations (C.6) and f6.s; 

(i,l) 


c - of particular interest. These equations governing R \ ’ / i 

( ■ 2 ' A ^ 

’ ' (x,x) which reduce for x = x to the Reynolds stress anc 


(x,x) and 
d the turbu- 
lent heat flux, respectively, quantitatively describe their nature as 
mechanisms driving the shear stress and the heat flow at the molecular level, 
in principle, these equations allow us to determine evolution of the 
Reynolds stress and the turbulent heat transfer beginning with fluctuations 
at the molecular lever and ending with fully developed turbulence. 


7. MACROSCOPIC DERIVATION OF 
CORRELATION EQUATIONS 

A we have shown in Section 5, the one-particle equations for a com- 
pressible tu*"!” !.e.nc flow which were originally derived on the basis of 
kinetic theory, could be derived also from the phenomenological equations 
of fluid dynamiics. Tl. are fore, it might well be expected that our two-point 
equations, viz., equations (6.4) through (6.9), (6.8) , and (6.19) , should 

similarly be interpreted in terms of fluid-dynamic equations (5.1) through 

o o o 

(5.3) governing "fluctuating" quantities p, u, and p. We note chat these 
equations are c dered as deterministic only in the sense that they have 
as many different solutions as number of different choices of the initial 
values I , (s: 1 , . . ,N) at the molecular level. Therefore it 

would be more realistic to categorize the equations as stochastic because 
specifying 6N inicial values is not practicable. Then, transition to the 
deterministic scheme of gasdynamics is effected through averaging ac various 
hierarchy levels, the lowest of which has led to the one-particle equations 
(5.6) and (5.9). In order to construct the two-particle hierarchy out of 
the stochastic fluid-dynamic equations we proceed as follows: let all the 

fluid-dynamic correlations be decomposed in the form 
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This expression being substituted in equation (6.6) we have, except for the 
term.s of thermal agitations. 
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First we note, in view of equations (7.1), (3.2) and (1.15) , that 
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we have, then, from (3.12), 


= Aa 
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With these linear relations for fluctuations, equation (7.2) is w-ritten in 
an equivalent form 



(7.5) 


where L. has been defined by (5.2). By a similar procedure we can show 
that the whole system of correlation equations (6.4) through (6.9) is simplv 
equivalent to 
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A set of equations = 0 and = 0, consisting of 5 and 15 indepen- 

dent equations, respectively, forms a complete set to determine 20 unknowns, 

(*J K ) 

namely, p , m^ or w^ , p, and the lovvest 15 components of the tensor R ’ 

For incompressible flows the equations with component - 4 (the energy 
equation) are deleted. On the other hand, p( = const.) and = 0) 

are not considered as dependent variables; accordingly we have 4 + 10 

4 

unknowns and the same number of equations. Tf, further, the flow is homo- 
geneous and isotropic, the first set of equations turns out to be trivial, 

whereas the second set is equivalent to the Karman-Howarth equation 

14 

describing the final stages of decay. 

Summarizing, the closure problem of the compressible, shear turbulence 
equations reduces, under the assumption of ternary molecular chaos, to 
solving equations (7.6) in conjunction with (4.3) through (4.5) supplemented 
by (4.12) and (4.13). These two systems are coupled through only two 
variables and which intervene in the latter groups of 

equations. Despite their seeming complexity, the closure equaticns may be 
solved by the method of separation of variables in (x,t) and (x,t) as is 
inferred from their rearranged form (7.6). The practicability of this 

4 

method has been ensured in the case of incompressible flows. Obtaining 
solutions, however, is beyond the scope of the present report. 

Nielst.n Engineering & Research, Inc. 

Mountain View, California 
February 1976 
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TABLE I 


PHYSICAL IMPLICATIONS OF ^OR (J,K) < (2,z) 




Expansion 

Coefficients 

o(O.O) 




„( 2 . 0 ) 




Incompressible Physical Meaning 

Flows I 

0 Density-density correlation 

n Turbulent correction to the mean 

velocity as distinguished from 
the mean mass-flux density divide 
hy the density (cf. eq. (35)) 

Nonvanishing Reynolds stress (of. eq. (54)) 

0 Turbulent correction to . the mean 

° temperature as distinguished from 

the pressure divided by the density 
(cf . eq. (36) ) 

Nonvanishing* Tv.rbulent heat flux density (of. 

eq. (37) ) 


(3PR)-=r(2>^) Nonvanishing* Temperature-tempera ture correlation 

*But has no practioal meaning (aeoause the temperature is a redundant 
variable in incompressible flow equations. 


